Odd-Integer Quantum Hall Effect in Graphene: Interaction and Disorder Effects 
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We study the competition between the long-range Coulomb interaction, disorder scattering, and 
lattice effects in the integer quantum Hall effect (IQHE) in graphene. By direct transport calcu- 
lations, both V = 1 and = 3 IQHE states are revealed in the lowest two Dirac Landau levels. 
However, the critical disorder strength above which the = 3 IQHE is destroyed is much smaller 
than that for the u = 1 IQHE, which may explain the absence of a = 3 plateau in recent exper- 
iments. While the excitation spectrum in the IQHE phase is gapless within numerical finite-size 
analysis, we do find and determine a mobility gap, which characterizes the energy scale of the sta- 
bility of the IQHE. Furthermore, we demonstrate that the v = 1 IQHE state is a Dirac valley and 
sublattice polarized Ising pseudospin ferromagnet, while the v = 3 state is an xy plane polarized 
pseudospin ferromagnet. 
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A number of dramatic recent experiments [3, [3, S, I3| 
have demonstrated the Dirac-like character of the low- 
energy electrons in graphene, a single monolayer film of 
carbon exfohated from graphite. In a relatively weak 
magnetic field, where the Zeeman splitting is negligible, 
an unconventional quantization of the Hall conductivity 
is observed, axy = i^j;- with = 4(fc -I- i) and k an in- 
teger 0, 1^. This can be ascribed to the Berry phase 
anomaly at the Dirac points 0, 0, d, 0,0, S] and the 
four-fold spin and sublattice symmetry [9| (pseudospin) 
degeneracies of the Landau levels (LLs). Interestingly, 
additional odd-integer = ±1 Hall plateaus together 
with even-integer — ±2, ±4... Hall plateaus were ob- 
served in a recent experiment [lo| by using a strong mag- 
netic field. A magnetic field which is sufhciently strong 
to lift the spin degeneracy of the LLs is expected to pro- 
duce the quantization rule f = 2k, as illustrated in Fig. 
1, which explains only the even-integer Hall plateaus. 

The even parity of v is assured in the clean, non- 
interacting limit by the valley degeneracy of the two 
Dirac points, which in turn is protected by the point- 
group symmetry of ideal graphene. The odd-integer 
quantum Hall effect (IQHE) is considered by most au- 
thors to be caused by electron-electron interactions 111 . 

|ll|. These works obtain a pseudospin fer- 
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romagnetic (PFM) = 1 state 
associated with Haldane's repulsive pseudopotential [17| . 
based on the low-energy continuum two-valley Dirac 
fermion description. In the continuum limit, the point- 
group and spin-rotation symmetries of the material are 
elevated to a full SU(4) symmetry, which reduces to an 
SU(2) symmetry when Zeeman splitting is introduced. 
Using the Stoner criterion ll| , Nomura and MacDonald 
have obtained a phase diagram, where the v = 1 IQHE 
state has a much lower critical magnetic field than the 




FIG. 1: (color figure online) (a) A rectangular sample of 
graphene of size Lx x Ly. (b) Illustration of the electron 
density of states of the LLs (filled area) and the even-integer 
quantized plateaus of the Hall conductivity a^y (thick line) 
in the absence of electron interaction. 



v — 3 state for a given sample mobility. However, direc- 
tion of the SU(2) symmetry breaking (orientation of the 
PFM magnetization) is not determined from the contin- 
uum theory. It depends instead upon residual effects of 
the lattice, as addressed by Alicea and Fisher who 
obtained an easy-axis orientation corresponding to sub- 
lattice (charge density wave) order in the v = 1 state. 
Moreover, the energy gap measured in transport is also 
sensitive to disorder at the lattice scale. This is especially 
important here, because the low-energy excitations of the 
J/ = ±1 IQHE states may be gapless [13[, which may lead 
to a non-trivial energy scale characterizing the stability 
of the IQHE. When the higher odd-integer Hall plateaus 
with \v\ > 1 are observable is still controversial. To re- 
solve these issues, an exact account of the competition 
between the long-range Coulomb interaction, disorder, 
and lattice effects is desirable, but so far lacking. 
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In this Letter, we carry out exact diagonaUzation 
calculations in a honeycomb lattice model, which cap- 
tures all these effects naturally. Through direct trans- 
port calculations, we provide numerical evidence that the 
Coulomb interaction can induce the v = ±1 and ±3 Hall 
plateaus. It is shown that, when the disorder is relatively 
weak, a number of low-energy many-particle states carry 
a same constant Chern number, forming a mobility gap, 
which protects the IQHE. The critical disorder strength 
for the V = \ state, determined as the point where the 
mobility gap vanishes, is much greater than that for the 
V = i state, suggesting that the v = 3 IQHE may be 
observed experimentally if disorder scattering can be fur- 
ther suppressed. The i/ — 1 state is clearly demonstrated 
to be a pseudospin ferromagnet with Ising anisotropy in 
the weak disorder regime. Moreover, our energy spec- 
trum analysis indicates that a PFM order exists in the 
h' = 3 state with the easy axis polarized in_the xy plane, 
consistent with the theoretical suggestion 
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Our model Hamiltonian in a perpendicular field B is 
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where is the non-interacting Hamiltonian 
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FIG. 2: (color figure online) (a) The energy spectrum in the 
units Vao/i as a function of total pseudospin 25*2 calculated 
for a sample size of 96 x 96 and M = 4 x 96 at = 1 (iVe = 12) 
in the clean limit [W = 0). Here, the static Coulomb inter- 
action energy has been included, and the irrelevant constant 
Zeeman energy shift has been omitted, (b) Anisotropic gap 
energy A a as a function of ao/l, where the sample sizes range 
from 36 x 36 to 120 x 120 for iVe = 12 and A^e = 9. The dotted 
line represents a parabolic fit to the data, (c) Excitation gap 
Ae as a function of 1/Ne with the dotted line as a linear fit. 
The values of M are chosen so that v = 1 or Ns = N„. 
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and the second term in Eq.(IT]) is the Coulomb interaction. 
Here, ni = c^ia'^ia the electron number operator on 
site i, tij = te'°'^ is the electron hopping amplitude be- 
tween neighboring sites in the presence of a magnetic flux 
(j) = J2o^v ~ W P^"- hexagon with M an integer, 
gaB is the Zeeman coupling energy with cr = ±1 for 
electron spin parallel and antiparallel to B, and Wi is a 
random on-site potential uniformly distributed between 
[~W/2, W/2], accounting for nonmagnetic disorder. De- 
noting the nearest neighbor carbon-carbon distance by 
flo, the magnetic length £ defined as usual is given by 

We first diagonalize the noninteracting Hamiltonian 
Hq on a rectangular sample (Fig. la), and obtain the 
complete set of single-particle wave functions of Hq. For 
the range of fields and disorder strengths considered here, 
the LL broadening from disorder scattering is always 
small compared to the LL spacing, and so the states asso- 
ciated to a given LL are clearly identifiable. We assume 
that the magnetic field is strong enough to cause com- 
plete splitting of the LLs for two spin directions. The 
total degeneracy of each LL near band center is denoted 
as 2Ns for each spin, i.e., Ns{= ^m" ) is the degeneracy 
for each Dirac component. We define Nf. as the electron 
number in the highest occupied LL - the rfi^ - such that 
the number of electrons counted from the band center 
is 2nNs + Ne, with < < 2Ns. The filling num- 



ber is = 2?! -|- Ne/Ng. Because of full spin polariza- 
tion, the relevant matrix elements of the Coulomb in- 
teraction are those with i ^ j, which are taken to be 
;7(R, - Rj) = Vao/\Ri - Rj|. The Coulomb interac- 
tion is projected into the n-th LL, and the many-particle 
wavefunctions are solved exactly in the subspace of the 
LL. 

For filling number < ly < 2, the Fermi energy is lo- 
cated inside the lowest n — LL. Denoting by A, B the 
two sublattices of sites, the z-component of the pseu- 



dospin Sz is expressed as 2Sz 

(in h = 1 units), which is conserved as the central LL 
eigenstates can be chosen to have support only on one of 
the two sublattices (the correction from lattice model is 
smaller than 10~® for system sizes that we consider). In 
Fig. 2a, we show the calculated many-particle low-energy 
spectrum at = 1 for = as a function of 2Sz, where 
Lx = Ly = 96, and M = 4 x 96. Periodic boundary 
conditions are imposed in the x and y-dircctions. 

In Fig. 2a, the lowest row of A'e -I- 1 energies corre- 
sponds to PFM states for N^ + l different eigenvalues of 
2Sz between -Ne and Ne. The two with 2Sz = Ne and 
—Ne have the lowest-energy, with intermediate values 
—Ne < 2Sz < Ne exhibiting higher energies. Clearly, this 
result suggests the presence of pseudospin anisotropy, 
with the z— axis as the easy axis 12|. In more physical 
terms, the favored 2Sz = iNe values represent charge 
ordered states with electrons occupying only one sublat- 
tice. We can define an anisotropic energy A a equal to 
the energy difference between the lowest eigenenergies at 
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FIG. 3; (color figure online) (a) The energy spectrum in the 
units Vat^ll as a function of total momentum q in Landau 
gauge calculated for a sample size 96 x 96 and M = 4 x 96 at 
= 3 (A^e = iVs) in the clean limit (VK = 0), (b) the (A^e + 1) 
lowest eigenenergies, as indicated by circle in (a), fitted by a 
parabolic function of total pseudospin 2Sz, and (c) ground- 
state energy for A'^e = 10 and = 12 calculated for different 
sample sizes and M with the dotted lines as a guide to the eye. 
In (c), for A^e = 10, the cross and diamond symbols indicate 
Ising and xy plane PFM states, respectively, and all systems 
with A^e = 12 (circles) are in the xy plane PFM state. 



25^ = -A^e and at 2S:, = -(TV^ - 2). Aa calculated 
for several different sample sizes is shown in Fig. 2b as 
a function of a^jl. The data can be well fitted by a 
parabolic function cx: [a^jt)^, which vanishes in the 
continuum limit faster than the characteristic Coulomb 
energy Vao/£. This is consistent with the interpretation 
of the pseudospin anisotropy as arising from corrections 
due to lattice effects, resulting in an additional ao/£ sup- 
pression factor. 

In Fig. 2a, we also see a small energy gap Ae between 
the PFM ground state and the lowest excited state in the 
second lowest row. We calculated Ae for different values 
of electron number A'e from A"e = 6 up to 24, as plotted 
in Fig. 2c as a function oil/Ne, where the magnetic flux 
strength 1/M is chosen to be nearly constant at different 
A^e, such that changes proportionally with the sample 
size X Ly. The data can be roughly fitted by a linear 
relation Ae oc l/N^- We note that in the absence of 
anisotropy, such gapless Ae ^ l/L"^ ^ behavior 
would be expected for the first excited pseudospin-wave 
states with \q\ ^ 1/L. Though the Ising anisotropy would 
be expected to introduce a gap, the observed behavior 
is probably consistent with the rather small anisotropy 
energy (note the scale in Fig. 2b). 

We have also carried out a spectral analysis for filling 
number v = where half states in the n = 1 LL are filled. 
Though in the continuum limit, the absence of coupling 
between valleys means that the pseudospin is conserved 
in this LL, there is no obvious S, conservation on the lat- 
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FIG. 4: (color figure online) (a) Hall conductivity a^y for 
disorder strength W = O.lt and three different interaction 
strengths V, averaged over 40 disorder configurations. Here, 
Lx = Ly — 54, M = 3 X 54, and the error bars stand for 
the standard deviation due to disorder average, (b) and (c) 
are the probability distributions of the ground-state Chern 
number for V = 0.25t and V = 0, respectively. 



tice analogous to the n = case. We show in Fig. 3a the 
low-energy spectrum in each total momentum q sector 
for pure system W — and system size L^ = Ly — 96. 
Interestingly, the lowest + 1 energies are all in the 
q = Ne/2 (in units of ^T^) sector with no double oc- 
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cupancy of any of the pseudospin doublets. Thus they 
are low-energy spin excitations, which can be fitted into 
I\E = {En - Eo) = aSl (with a > 0) as shown in Fig. 
3b. This suggests that the nondegenerate ground state 
has Sz = 0, and is an xy plane polarized PFM state, with 
strong valley mixing. We have further checked a number 
of system sizes between 24 x 24 to 200 x 200, and found 
that the xy plane polarized state is always the ground 
state as long as both Lx and Ly are commensurate with 3 
(that includes all the systems with Ne = 12). Otherwise, 
an Ising PFM state is found to be favorable, as shown 
in Fig. 3c. This strong systematic finite-size effect can 
be understood from the graphene band structure, since 
valley mixing implies order at the wavevector connecting 
the two Dirac points, and hence period 3 modulations in 
both lattice directions [l^. Indeed E^) shows an oscilla- 
tion with an upturn at Ising points, indicating frustration 
of the modulations in the energetically preferred xy PFM 
state. The xy plane PFM state is expected to become 
the ground state for = 3 at the thermodynamic limit. 
The charge density is uniform in the xy plane state with 
vanishing charge current on each lattice bond. Interest- 
ingly, in the Ising state, we observe lattice-scale charge 
currents circulating around one third of the hexagons in 
the pattern predicted by Alicea and Fisher [l^ . 

Given that any gap for the v = 1 IQHE is small enough 
to be numerically unresolvable, it is important to di- 
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FIG. 5: (color figure online) (a) Calculated Chern numbers of 
60 low-energy eigenstates as a function of E„ — Eq with En 
the n-th eigenenergy at u — 1, for V = 0.5t, W — 0.8V and 
10 random disorder configurations, (b) Critical energy Ec for 
filling numbers v = 1 (squares) and v = 3 (triangles) as func- 
tions of normalized disorder strength W/V, where the error 
bars are the mean deviation of Ec due to disorder sampling. 



diminishes to zero, where the v = 1 IQHE is destroyed. 

By similar calculations, we find that odd IQHE can 
also occur in higher LLs, in consistence with the xy plane 
PFM order. The calculated phase diagram ioi ly — 3 
IQHE in the n — 1 LL is shown in Fig. 5b (triangles). 
The V ^ 3 IQHE is less stable than the 1 IQHE, 
with a critical disorder strength Wc — 0-QV about one 
third of that for v = 1. This may explain the observation 
of the u = 1 but not v ~ 3 plateau in experiment [l3| ■ 
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rectly demonstrate its robustness to disorder. We now 
calculate the Hall conductivity axy, which can be ex- 
pressed in terms of the ensemble average of the Chern 
number 2^ 21| Co of the ground state as axy = ^(Co). 
In Fig. 4a, the calculated a^y, averaged over 40 random 
disorder configurations, is shown as a function of filling 
number for a weak disorder strength W = O.lt. In the ab- 
sence of Coulomb interaction (F = 0), axy increases con- 
tinuously with v, without showing a quantized plateau 
around ^ I. However, as the interaction is switched 
on, a quantized Hall plateau appears around 1^=1. In 
Fig. 4b, the Chern number distribution for V — 0.25t at 
filling numbers = 0, i, ■ • ■ 2 is shown. Near integer fill- 
ing numbers 1 and 2, the Chern number takes constant 
values Co = 1 and Co = 2 for all disorder configurations 
without fluctuations, corresponding to the i> — 1 and 
1^ — 2 IQHE plateaus in Fig. 4a, respectively. For V = 0, 
as shown in Fig. 4c, various Chern numbers, Co=0, 1 
and 2, merge together in the middle region, resulting in 
a plateau-metal transition. 

We now study the thermal stability of the odd IQHE 
by also considering the excited states. In Fig. 5a, we show 
the Chern numbers of 60 lowest eigenstates calculated at 

~ 1 ior Lx — Ly — 60 and iVe = 10 as a function of 
En — Eq. The Chern numbers for 10 random disorder 
configurations of strength W = 0.8V are represented by 
different symbols. We see that the Chern numbers of low- 
energy eigenstates with £'„ — Eq smaller than a critical 
energy Ec always take a constant value C„ — 1, indicat- 
ing localization for these states and a mobility gap (which 
is directly related to the activation gap) of order Ec [IH . 
The calculated Ec as a function of W/V for V = 0.5t is 
shown in Fig. 5b (squares). For W > Wc — 2.0V, Ec 
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